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Abstract 

A review of gravitational and hydrodynamical processes during forma- 
tion of clusters and evolution of galaxies is given. Early, at the advent 
of N-body computer simulations, the importance of tidal fields in galaxy 
encounters has been recognized. Orbits are crowded due to tides along 
spiral arms, where the star formation is enhanced. Low relative velocity 
encounters lead to galaxy mergers. The central dominating galaxies in 
future clusters form before the clusters in a merging process in galaxy 
groups. Galaxy clusters are composed in a hierarchical scenario due to 
relaxation processes between galaxies and galaxy groups. As soon as the 
overall cluster gravitational potential is built, high speed galaxy versus 
galaxy encounters start to play a role. These harassment events grad- 
ually thicken and shorten spiral galaxy disks leading to the formation 
of SO galaxies and ellipticals. Another aspect of the high speed motion 
in the hot and diluted Intracluster Medium (ICM) is the ram pressure 
exerted on the Interstellar Matter (ISM) leading to stripping of the ISM 
from parent spirals. The combinations of tides and ram pressure strip- 
ping efficiently removes the gas from spirals, quenching the star forma- 
tion in galactic disks, while triggering it in the tidal arms and at the 
leading edge of gaseous disk. Gas stripping from disks transports the 
metals to the ICM. In some cases, the gas extracted from the galactic 
disks becomes self- gravitating forming tidal dwarf galaxies. 

Star formation (SF) is another important ingredient in the evolu- 
tion of galaxies. Young stars provide the energy, mass feedback and the 
metals to the ISM. SF also drives the supersonic turbulence and trig- 
gers SF at other places in the galaxy. Examples of supershells in the 
ISM, resultant from the evolution of stellar cluster, are given. Struc- 
tures, supershells, filaments and sheets are also produced when the ISM 
is compressed during galaxy collisions and when high velocity clouds 
(HVC) fall and ram through gaseous galaxy disks. In some cases, the 
compressed structures become self- gravitating and fragment to produce 
clumps. When a low density medium is compressed, the clump masses 
and sizes correspond to those of giant molecular clouds (GMC). When 
already existing CMC's are compressed they are driven to collapse form- 
ing massive super-star clusters. The fragmentation process in the dense 
environment of GMC's produces fragments of stellar masses with the 
mass function having a slope similar to the stellar Initial Mass Func- 
tion. 
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1 Introduction 



The internal evolution of isolated galaxies is the subject of many studies. The 
evolution of stellar disks depends on the ratio between the radial component 
of the velocity dispersion or multiplied with the local epicyclic frequency k, 
which is an analog of pressure in the stellar disk, and the disk surface density 
S multiplied with the constant of gravity G, which gives the local gravitational 
force in the flat disk. The above quantities are combined in the Toomre (1964, 
1977) Q parameter 

Q= aRK (I) 

to evaluate the local stability of a rotating, self-gravitating disks: for Q < 1 
the disk is locally unstable and forms large scale deviations from the axial 
symmetry - spiral arms. A purely stellar disk heats up, when individual stars 
scatter on spiral arms, and self-stabilizes. 

If the fraction of the disk kinetic energy in random motions is small enough, 
the bar forms in its central part (Ostriker & Peebles, 1973). The bar exchanges 
the angular momentum with stars, which results in radial redistribution of 
mass leading to the formation of a central mass concentration. The growing 
central mass concentration partially dissolves the bar itself (Shen & Sellwood, 
2004). 

A dissipative gaseous interstellar medium (ISM) is an additional compo- 
nent of the disk moving in a common gravitational field with stars. Energy 
dissipation, supersonic shocks and collisions of interstellar clouds reduce ran- 
dom motions of the ISM agitating spiral-like instabilities, which gradually 
heat the stellar disk. In an isolated spiral galaxy, the heating by a bar and 
spiral arms increases the stellar velocity dispersion by less than a factor of 2 
(Bournaud et al., 2004). 

The ISM complements the influence of the central mass concentration with 
viscous forces that shift systematically the gas flow-lines relative to the stellar 
orbits. The evolution of a bar and spiral structure in an isolated galaxy, their 
growth and dissolution, is driven by a combination of the increase in stellar ve- 
locity dispersion due to scattering of individual stars on large scale deviations 
from the axial symmetry, by the growth of the central mass concentration 
and by the torque between the gaseous and stellar components (Bournaud & 
Combes, 2004). 

Star formation, the mass recycling between stars and ISM, and the energy 
feedback are also recognized as fundamental processes. They are included in 
closed box models the evolution of isolated galaxies (Jungwiert et al. 2001, 
2004). Several issues are addressed, remaining as open question if, when, and 
how much the box has to be open to gas infall, environmental effects, and 
galaxy major and minor mergers: 

• Persistence of star formation. With the present rate, star formation 
would consume all the gas within less than a few Gyr. Fresh gas supplies 
are needed to keep the star formation active. 
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• The metalicity distribution and chemical evolution of the disk. The 
G-dwarf problem requires gas infall. 

• The existence of a thin gaseous disk, where recurrently the instabilities 
operate also needs accreted fresh gas. 

• Thin disks are destroyed in minor merger and harrasment events. They 
have to be recreated. 

• Renewal of bars driving the mass to the galactic center also requires gas 
accretion. 

In this review, the tides, major, intermediate and minor mergers events, 
harrasment events, gas stripping, star formation and feedback are described. 
We discuss the relevance of these processes in the formation, and evolution of 
galaxies, and in the evolution of galaxy groups and clusters. 

2 Tides 

The importance of tides has been shown in early numerical N-body simula- 
tions by Toomre & Toomre (1972). Their restricted N-body simulations of an 
encounter between two spiral galaxies led to the basic features. As the galax- 
ies pass by each other, tidal forces provide the stars and ISM of the disks 
with sufficient energy to escape the inner potential well. Two streams of disk 
material appear on both near and far side of the disks relative to the position 
of the peri-center. In the near side, stars and gas form a tidal bridge between 
the disks while the far side material forms long tidal tails. The length and 
the shape of the tails is a sensitive function of the relative velocity and of the 
geometry of the encounter: the most effective are the low relative velocity pro- 
grade encounters, where the spin and angular momentum vectors are aligned. 
In this case, the momentum and energy exchange in resonant orbits is most 
effective, leading to long tidal arms and a close interconnection between the 
central parts of the original disks, to merge in the future. 

The evolution of tidal debris is described by Mihos (2004). Most of the 
remnant material remains bound. It follows elliptical orbits and only a small 
fraction of the orbits are unbound. The fraction of gravitationally bound 
material is proportional to the extend of the dark matter halo: large halos 
result in less unbound material. The bound tidal debris return after some time 
to the galaxy disk in the form of gas infall and high velocity clouds, which 
may trigger the star formation. A fraction of the disk gas on unbound orbits 
becomes part of the Intracluster medium (ICM). The metalicity of the ICM, 
which is about 1 /4 Z & , gives an evidence of the connection between galaxy 
disks and ICM. The slow encounters are the most effective in driving the gas 
of the original disks inward to the center of the future galaxy produced by the 
merger, which triggers a nuclear starburst. 
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3 Major mergers 



Computer simulations of low velocity encounters of spiral galaxies with similar 
masses - major merger events - produce remnants that are in surprisingly good 
agreement with the observed shapes, density profiles and velocity distribution 
of the observed giant elliptical galaxies. The mass ratio of the progenitor disks 
determines the global properties of the remnant (Burkert & Naab, 2003a, 
b). Giant elliptical galaxies can be subdivided into disky and boxy types 
(Naab & Burkert, 2003): the less massive disky giant ellipticals show disk- 
to-bulge rations of SO galaxies and rotation in their outer parts. They can 
be understood as resulting from major mergers with mass ratios from 1/1 
to 1/4. The more massive boxy giant ellipticals form by 1 : 1 mergers with 
special initial orientations only, which weakens the disk merger scenario as 
the possible formation mechanism. Other processes have to be added, such 
as the star formation and the energy dissipation in the gaseous components 
of progenitor disks. 

The present day giant elliptical galaxies may have been formed by major 
mergers of gas- rich galaxies with a subsequent starburst, or by mergers of 
less gas-rich galaxies, which is more common in galaxy clusters. The second 
scenario assumes that a fraction of stars forms before ellipticals in smaller gas- 
rich galaxies. In the ACDM cosmological simulations the number of mergers 
varies with z. The giant haloes of future elliptical exist already at z = 6 
with the relatively small variation of mass at z < 6. Giant ellipticals exist at 
redshifts z = 3 with about 50% of stellar mass, however, a typical cD galaxy 
has suffered significant merging events even at redshifts z < 1 (Gao et al., 
2004). Bimodal distribution of metalicity observed in elliptical galaxies is a 
product of gas-rich mergers, when the globular clusters form in the tidal arms 
(Li at al. 2004). 

The Antennae, the colliding pair of galaxies NGC 4038/39, provide an 
nearby example of an early phase of a low velocity encounter between galaxies. 
The velocity field of Antennae galaxies, showing the details of the interaction, 
has been measured by Amram et al. (1992). Later, when the collision partners 
will not be distinguished any more, it should result in a field elliptical galaxy. 
Another example of a slightly more advanced ongoing merger event, is the 
Hickson compact group 31 (Amram et al. 2004). HCG 31 is a group in early 
phase of merger growing through slow and continuous acquisition of galaxies 
from the associated environment. One of the best examples of remnants after 
a recent merger event is the giant elliptical galaxy NGC 5018. Its inner part 
shows a uniform ~ 3 Gyr old stellar population presumably produced in a 
merger induced starburst (Buson et al., 2004). 

4 Intermediate and minor mergers 

Collisions of intermediate (1/5 - 1/10) mass ratio partners show the inter- 
mediate merger events, which are explored by Bournaud et al. (2004). They 
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result in a hybrid system with spiral-like morphology and elliptical-like galaxy 
kinematics similar to some of the objects identified by Chitre & Jog (2002). 

Minor mergers, i.e. merger of a galaxy with a satellite 1/10 or less massive 
than the galaxy itself, are discussed in the context of formation of thick disk 
and of disk globular clusters (Bekki & Chiba, 2002, Bertschik & Burkert, 
2003) . The disk globular clusters are formed in the high-pressure dense central 
region of the gas-rich dwarf galaxy, which is compressed in a tidal interaction 
with the thin disk of the spiral galaxy, and later they are dispersed in the disk 
region, when the original dwarf merges with the galaxy. The origin of thick 
galactic disks may be the result of the same minor merger event producing 
the disk globular clusters. 

5 Evolution of galaxies in clusters 

An overabundance of spiral galaxies and an under abundance of SO galaxies in 
high-redshift clusters, when it is compared to low z clusters, is a consequences 
of galaxy merger events and of other environmental effects. In ADCM cos- 
mological models with Oo = 0.4 and fl\ = 0.6, the galaxy interactions and 
tidal effects play a role mainly at intermediate redshifts 0.5 < z < 5 (Gnedin, 

2003) . Galaxies entering clusters with low relative velocities merge their halos 
with the cluster and their subsequent dynamical evolution is due to perturba- 
tions along their orbit. The infalling galaxy groups experience merger events 
of their members. After virialization, when the galaxy velocity in the clus- 
ter increases to a few 10 3 km s _1 the mergers are suppressed and high-speed 
galaxy encounters and interactions with the intracluster medium (ICM) start 
play a role. 

The ICM consists of hot (~ 10 7 K) diluted (10~ 3 - 10~ 4 cm" 3 ) gas, which 
is detected in X-ray observations. The optical/infrared observations (Ren- 
zini, 1997) monitor the stellar component of galaxies, and also their chemical 
composition. The bulk of the cluster light is produced in bright giant ellipti- 
cal galaxies and in galactic bulges by intermediate and low mass stars. The 
content of Fe in clusters scales with their total light and the abundances of 
various elements does not seems to change from cluster to cluster (Renzini, 

2004) . If the metals are mainly produced in SN Type la and SN Type II, the 
constancy of metallic ratios implies the same universal ration of the two types 
of SN. It may be interpreted as a sign that majority of cluster metals have 
been produced in the dominating giant galaxies, and that the stars formation 
holds a universal initial mass function (IMF). It appears that the total mass 
of metals in the ICM supersedes the total mass of metals in stars (Renzini, 
2004). The galaxies lose more metals to the ICM than they are able to retain. 
Their total abundance does not correspond to the actual SN rates, the obser- 
vations could be reproduced if the SN Type la in ellipticals have 5-10 time 
larger rates in the past compared to present and if the slope of the IMF is not 
too shallow: Salpeter (1955) IMF with the power law slope -2.35 matches the 
requirements. 



5 



How the metals were transported to the ICM? Potential mechanisms are 
the ICM ram pressure stripping of galactic ISMs, which may become more ef- 
fective in combination with gravitational tides, or the star formation feedback 
connected to early stellar winds driven by the starburst forming majority of 
the stellar galaxy itself. An evidence for early stellar winds serve the Lyman- 
break galaxies (Pettini et al., 2003). Another possibility are late local winds 
due to massive starbursts (Heckmann, 2003) or flows driven by a declining 
rate of SN type la (Ciotti et al., 1991). Both, gas extraction by stripping and 
gas ejection by winds will be described below. 

6 Galaxy harrasment 

Moore et al. (1996) discuss the heating effects influencing sizes and profiles of 
dark matter haloes of galaxies in clusters: tidal heating of haloes on eccentric 
orbits and impulsive heating from rapid encounters between haloes. N-body 
simulations show that both processes restrict the halo sizes. Frequent high 
speed galaxy encounters in clusters, galaxy harassment, drive the morpho- 
logical transformation of galaxies in clusters. Moore et al. (1995) show the 
dramatic evolution, which happened in clusters during the recent ~ 5 Gyr. 
Both, young clusters at z ~ 0.4 and old nearby clusters have central dominant 
elliptical galaxies, which have been formed before clusters. The young clusters 
are populated by many spiral galaxies, but the old, large virialised systems, 
have inside of a deeper potential valley all other galaxy types, like SO, dwarf 
Ellipticals, dwarf SO, dwarf Spheroidals and Ultra Compact Dwarfs (Moore, 
2003). Numerical simulations show that if the cluster velocity dispersion is 
more than a few times the internal velocity dispersion of galaxies, they do not 
merge in encounters. The observed morphological change is due to impulsive 
interactions during high speed encounters restricting the dark matter haloes 
and pumping the kinetic energy into disks of spirals changing them to SO, dE, 
dSO, sSph, UCD galaxies. Hydrodynamical processes, like gas stripping, also 
matters, and help with the change, particularly in relation to star-formation 
and ICM. 

7 Gas stripping 

Gunn & Gott (1972) in a description of events during and after the collapse, 
or formation of a galaxy cluster, predict the formation of hot and diluted intr- 
acluster medium (ICM) produced in shock randomization of infalling gaseous 
debris. The interstellar medium (ISM) in a galaxy moving through the ICM 
feels the ram pressure P r , which is proportional to 

Pr ~ Picmv 2 , (2) 

where pic m is the density of the ICM and v is the velocity of the galaxy rela- 
tive to the ICM. The ISM may be stripped away the parent galaxy if at given 
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distance from the galactic center the ram pressure exceeds the gravitational 
restoring force F r 

F r = 27rGS s E/5M, (3) 

where E s is the surface density of stars, Ejsm is the surface density of the 
ISM and G is the constant of gravity. 

The ram pressure gas stripping influences the presence of the ISM, which 
may affect the star formation rates in galaxies. It is still an open question, 
if the observation of Butcher & Oemler (1978, 1984), showing that distant 
clusters contain a far higher fraction of blue star forming galaxies than their 
near-by counter-parts, can be explained with the ICM ram pressure stripping 
of the ISM from galaxies. 

A mechanism that suppresses the star formation in local galaxy clus- 
ters is apparently connected to the morphological transformation of galaxies: 
SO galaxies are under-represented in distant galaxy clusters, luminous spiral 
galaxies are in deficit near the centers of local galaxy clusters (Dressier et al., 
1997; Couch et al. 1998). This suggests that the cluster environment removes 
the ISM from galaxies, suppressing star formation and transforming spirals 
to SO's. 

Numerical simulations of ram pressure gas stripping using 3-dimensional 
SPH/N-body code have been preformed by Abadi et al. (1999). They confirm 
the predictions of Gunn & Gott (1972) that the radius to which the gas is 
removed from the parent galaxy depends on the relation of the ram pressure 
to the restoring force. But in any case a substantial part of the cold gas 
remains sufficiently bound to the stellar disk. The star formation rate is 
reduced by a factor of 2 only, which brings them to a conclusion that the 
simple ram pressure stripping does not adequately explain the sharp decline 
of star formation seen in Butcher-Oemler effect. 

Volmer et al. (2001) show in simulations using sticky particles that the gas 
stripping is rather sensitive to the galaxy orbit, particularly to the minimum 
distance to the cluster center, and also to the orientation of the galactic disk 
relative to the orbit inside the cluster. In some cases the removed ISM is 
re-accreted and it falls back to the galactic disk, possibly triggering star for- 
mation in the central part of the disk within the remaining gas. Simulations 
of ram pressure gas stripping along an orbit in a cluster by Jachym (2004) give 
similar result: the orbital parameters and the orientation of the galaxy are 
important. The gas is more effectively removed when the pole of the galactic 
plane is near the direction of the galaxy orbital motion in a cluster (Fig. 
compared to the situation when the galaxy moves edge-on along its orbit. 

8 Gas stripping and tides 

Recently reviews on the efficiency of the stripping mechanism were given by 
van Gorkom (2004) and Combes (2004). The theoretical considerations stress 
the effects of viscosity and thermal conduction including Kelvin-Helmholtz 
instability (Nulsen, 1982). The turbulence, shells and supershells formed by 
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Figure 1: Simulations of ram pressure gas stripping (Jachym, 2004). Edge-on 
orientation of the galaxy plane relative to the orbit (left panel) is compared 
to the pole-on orientation of the galaxy plane (right panel) . 

star formation increase the effect of the ram pressure stripping due to erosion 
at the edges of ISM holes (Bureau & Carignan, 2002). Ram pressure also 
contributes to the overpressure inside of expanding structures increasing their 
sizes. 

A proto-type of overwhelming gas stripping in the Virgo cluster is the 
galaxy NGC 4522 with a truncated HI disk, enhanced star formation in the 
central region, extraplanar gas, while undisturbed stellar disk. The gas dis- 
tribution suggests an ICM-ISM interaction, while the undisturbed disk rules 
out a gravitational interaction (Kenney & Koopmann, 1999). 

The transformation of galaxies in dense environment results from the com- 
bined action of gas stripping and gravitational tides. C 153 galaxy in the 
cluster Abel 2125 shows an ongoing gas stripping: A tail of ionized gas is seen 
in [0 II] emission, which extends at least 70 kpc toward the cluster core along 
C 153 orbit, coincide with a soft X-ray feature seen in the Chandra obser- 
vations (Keel et al., 2004). At the same time the HST optical picture shows 
clumpy morphology, including luminous star-forming complexes and chaotic 
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dust features. The perturbed stellar disk with enhanced star formation ac- 
tivity suggests a possibility that a burst of star formation has been initiated 
during the close passage of C 153 to the cluster center. 

More examples of the stripped dwarf irregular galaxies (van Zee et al., 
2004) and of galaxies with truncated star formation disks in Virgo cluster 
(Koopmann & Kenney, 2004a, b) show both the gas stripping and gravita- 
tional interactions including the induced star formation, however, the dom- 
inant environmental effect on galaxies in clusters is the ram pressure gas 
stripping. 

9 Tidal Dwarf Galaxies 

The long tidal tails observed in many cases of interacting galaxies show mas- 
sive clumps of 1O 9 M (e.g. IC 1182, NGC 3561, NGC 4676, etc.). These 
massive blobs has been named Tidal Dwarf Galaxies (TDG), since they have 
galactic masses and the chemical composition corresponds to the recycled 
matter pushed out of the disks of the interacting partners. To become a long 
living independent galaxy they should be gravitationally bound systems (Due 
et al., 2000). The most prominent interacting system, the Antennae galaxies, 
has been studied by Mirabel et al. (1992), who describe a TDG candidate 
at the tip of their long tidal arm. Hibbard et al. (2001) also analyzed this 
concentration of gas and star forming regions, however, from observations it 
is very difficult to assess if it is gravitationally bound. 

To decide if a TDG candidate will be a new galaxy formed out of an 
interaction, several questions have to be addressed (Bournaud et al. 2004): 

• Are the blobs real concentrations in three dimensions? They may be 
just projection effects due to tidal arm geometry. 

• Are they kinematically decoupled from the tidal arms? Are they long 
living? 

• Do they contain dark matter? 

Bournaud et al. (2004) give the answer at least to the first question: The 
simulations of galaxy collisions provide shapes of tidal arms, which may be 
virtually observed from all the sides. They conclude that some observed 
10 9 Mq mass concentrations are real TDG candidates. Some of them are self- 
gravitating, but to decide on their future is still difficult with the current 
resolution of simulations. The star formation, energy and mass feedback have 
to be included in the future numerical experiments. The third problems on 
the content of the dark matter also remains open due to the uncertainty in 
the internal kinematics of TDG candidates. 

The simulations provide one more important conclusion (Due et al. 2004): 
the existence of 1O 9 M mass concentrations at the tips of tidal arms is rather 
sensitive to the extent and the density profile of the halo. It has to be ex- 
tended enough (~150 kpc from the center of a collision partner) so that the 
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collision happens within it. Then the flow lines of the perturbed gas from 
all galactocentric distances in the original disk concentrate at the tips of the 
tidal arms, giving a kinematical origin to the TDG. When the halo is too 
concentrated so that the collision happens at the distance, where the rotation 
curve already decreases, the perturbed gas populates all the tidal arm and 
there is no place where 10 9 M Q may gather. 

The kinematical gathering of stars and gas distinguishes TDG from super- 
star clusters (SSC), which are seen the interacting galaxies. SSC are not only 
less massive (10 5 — 10 7 M Q ) but they arise from gravitational instabilities in 
the stellar or gaseous components along the tidal arms. 

10 SMC, LMC and the Milky Way system 

The Small (SMC), and Large Magellanic clouds (LMC) and the Milky Way 
form the nearest interacting system of galaxies, where the gravitational and 
hydrodynamical processes can be studied. We see the result of a combination 
of gravitational tidal forces with gaseous ram pressure. The star formation 
and the mass and energy feedback is also involved as demonstrated with many 
expanding ISM shells and supershells. 

The high resolution HI surveys of the LMC (Kim et al. 1998, 1999), of 
the SMC (Staveley-Smith et al. 1997, Stanimirovich et al. 1999), of the 
Magellanic bridge region (Muller et al. 2003) and of all the system (Buns 
et al. 2004) have been preformed with the Australia Telescope Compact 
Array and with the Parkes radiotelescope. The following large-scale features 
are distinguished on the (1, b) integrated HI intensity maps: LMC, SMC, 
Magellanic Bridge joining the two clouds, Magellanic Stream starting at the 
SMC and following an almost polar plane passing less than about 10° of the 
galactic south pole and stretching more than 100°, the Leading Arms - the 
HI gas preceding the motion of the clouds, and the Interface Region - the HI 
between Magellanic Bridge and Magellanic Stream (Fig. [2}. All these features 
are also distinguished on the average radial velocity maps. The radial velocity 
changes smoothly from RVlsr = —400 km s^ 1 at the end of the Magellanic 
Bridge to +100 km s" 1 at the SMC, to +240 km s _1 at the LMC, while the 
Leading Arm does not show a clear gradient. Deprojection correcting for the 
solar motion reduces the velocity difference between the LMC and SMC from 
ARVlsr — 123 km s _1 to AKVqsr — 67 km s~ x . A further correction taking 
into account motions of the LMC relative to the MW reduces the difference 
to ARVlmcsr = 10 km s _1 . The HI average radial velocities show that 
the encounter between LMC and SMC happens at a small velocity not much 
larger that 10 km s" 1 , which makes the interaction rather long, giving the 
time to the gas to flow away from its parent cloud and form the observed 
features. 

The LMC HI disk seems to be compressed at the side opposite to the 
SMC, and the LMC shows rotation almost perpendicular to SMC - LMC 
direction. The LMC disk has a diameter of about 7.3 kpc with a rotation 
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curve rising rapidly to 55 km s _1 in the inner 1.5 kpc, more smoothly to a 
63 km s _1 peak at 2.4 kpc and declining thereafter (Kim et al., 1998). The 
SMC also shows rotation: SMC disk includes the bar-like feature about 4 kpc 
in extent, with a velocity gradient of about 100 km s _1 (Stanimirovic et al. 
2004). The rotation curve rises to about 60 km s^ 1 up to the turnover radius 
of 3 kpc. The velocity dispersion is high along the high column density axis 
of the Magellanic Bridge. LMC, Magellanic Bridge and SMC show similar 
velocity fields, which can result from the rotation of all the three partners 
with the same orientation along the axis LMC - SMC. The velocity field is 
rather broken at the southern part (in galactic coordinates) of the SMC and 
of the Magellanic Bridge. There, the Magellanic Stream and Interface Region 
start. At the other end, in front of the LMC, the Leading Arm can be split into 
three parts. Individual features in the Leading Arm show head-tail structure 
with the orientation along the direction of possible space motion of the LMC 
as it is given by Kroupa & Bastian (1997) and van der Marel (2001). 

The HI observation should be complemented with the studies of distribu- 
tion of planetary nebulae (Dopita et al. 1985), of carbon stars (Kunkel et 
al. 2000), and of Cepheids (Groenewegen 2000), and with the near- infrared 
star counts from the Two Micron All Sky Survey (2MASS) and the Deep 
Near- Infrared Southern Sky Survey (DENIS), van der Marel (2001) shows 
that LMC disk is not circular at larger radii, it is elongated in the direction 
of the Galactic centre, suggesting the influence of the tidal forces of the Milky 
Way. The data should be compared to models of the interaction in attempts 
not to be in contradiction. 

The N-body simulations modeling the gravitational interaction of the SMC 
with the LMC and Milky Way (Gardiner & Noguchi, 1996) show that the last 
two close encounters between the interaction partners, 1.5 and 0.2 Gyr ago, 
are able to explain many of the observed structures. Magellanic Stream and 
Leading Arm have been created as a consequence of the former close en- 
counter, 1.5 Gyr ago. The Magellanic Bridge and the Interface Region have 
been formed later during the last close encounter between LMC and SMC 
0.2 Gyr ago. The discovery of the Leading Arm stressed the importance of 
gravitational tides and questioned the role of ram pressure stripping (Putman 
et al. 1998). However, the separation between gaseous and stellar features, 
absence of stars in tidal Magellanic Bridge and Magellanic Stream, filamen- 
tary structures in the SMC and head-tail structures in the Leading Arm and 
in the Magellanic Stream show that also hydrodynamical forces like ram pres- 
sure gas stripping have to act there. Most of the gas in Magellanic Bridge 
and Leading Arm is coming from the SMC. The preencouter SMC, more that 
1.5 Gyr ago, must have a gas disk of about 10 kpc in diameter (Stanimirovic 
2004), which shrunk forming the Magellanic Stream and Leading Arm 1.5 
Gyr ago and Magellanic Bridge and Interface Region 0.2 Gyr ago. How- 
ever, still a substantial amount of angular momentum of the original disk is 
left corresponding to other simutions of dwarf galaxies (Mayer et al. 2001) 
demonstrating that the tidal stripping removes the angular momentum rather 
slowly, at the timescale of 10 Gyr. Consequently, the recent two encounters 
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GALACTIC Longitude 



Figure 2: LMC & SMC: an interacting system of dwarf galaxies (Briins et al. 
2004). 
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have not been able to remove a substantial part of the original angular mo- 
mentum making possible to see the rotation of the SMC at present times. The 
present rotation curve of the SMC shows that the dynamical mass within 4 
kpc is about 2.4 10 9 M Q three-quarters of which is stellar. The dark matter 
is not needed for an explanation of the rotation speeds in the SMC. 

Careful analysis of the HI distribution in the Magellanic Stream (Putman 
et al. 2003) suggests that the Magellanic Bridge is older that assumed above 
and the LMC and the SMC are bound together for at least two orbits. The 
dual filaments emanating from the SMC and from the Magellanic Bridge are 
of tidal origin and shaped by a small amount of ram pressure. 

Star formation and energy feedback from young stars are included in the 
so far most sophisticated N-body model of the LMC - SMC - Milky Way 
encounter by Yoshizawa & Noguchi (2003). This model agrees well with sev- 
eral observed features including the Magellanic Stream, which apparently has 
tidal origin, it reproduces the presence of young stars in the south-east wing 
of the SMC and it also reproduces the acceleration in the star formation ac- 
tivity, which is due to recent close encounters between the clouds. Some open 
and unsolved questions remain: bimodal or many peak distribution in the 
main gaseous body of the SMC and Magellanic Bridge remain to be inter- 
preted, it probably originates in the numerous expanding shells. The shells 
and supershells, if they trigger star formation, also remain an open question. 

The influence of other MW satellites is unclear. The tidal stream of HI 
clouds connected to disruption of the Sagittarius dwarf galaxy on the polar 
orbit aloud the Milky Way provides a possible explanation for the anomalous 
velocity distribution of HI clouds near the south galactic pole (Putman et al. 
2003) . The possibility of formation of the local group dwarf members including 
LMC and SMC out of the Milky Way encounter with the M31 galaxy is 
discussed (Sawa & Fujimoto, 2004). Another picture describes the interaction 
of the Fornax-LeoI-LeoII-Sculptor-Sextans stream with the Magellanic Stream 
causing the gas stripping from the Fornax (Dinescu et al. 2004). Kroupa et 
al. (2004) propose the origin of the whole Local Group with the local dwarf 
galaxies in a common great circle. 

11 Star formation, energy and mass feedback 

Star formation is a complex process of the gravitational collapse and frag- 
mentation, where the thermal and magnetic support competes with super- 
sonic shock waves and energy dissipation. The density increases by 20 orders 
of magnitude from that of a molecular clouds core. The interplay between 
gravity, magnetic forces, hydrodynamical processes, radiative transfer and 
chemistry happens in a turbulent interstellar medium. Supersonic flows form 
sheets and filaments involving mass concentrations, which in some cases are 
bound by self-gravity. Some of them collapse forming single or binary stars, 
others disperse. 

There are different sources powering the interstellar structures. On a small 
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scale the pre-main sequence stellar winds and stellar radiation, on somewhat 
larger scale the main sequence stellar winds, and on even larger scale the su- 
pernovae. Young stars pump energy back to the interstellar medium, which 
influences the conditions for further star formation. The energy released by 
young OB associations compresses the ambient medium into shells and super- 
shells, which may collapse and trigger new star formation. The energy feed- 
back triggering shell collapse and further star formation is a self-regulating 
mechanism of the galaxy evolution. 

11.1 The observation of shells 

Shells and supershells and holes in the HI distribution have been discovered in 
the Milky Way by Heiles (1979, 1984), in M31 by Brinks and Bajaja (1986), 
in M33 by Deul and Hartog (1990), in LMC by Kim et al. (1999), in SMC 
by Stanimirovic (1999), in HoII by Puche et at. (1992), in Ho I by Ott et 
al. (2001) and in IC 2574 by Walter and Brinks (1999). Most probably they 
are created by an energy release from massive stars, however, an alternative 
explanations, infall of high velocity clouds (Tenorio-Tagle and Bodenheimer, 
1988), or gamma ray bursts (Efremov et al., 1998; Loeb and Perma, 1998) 
has been invoked in some cases. The majority of the observed shells is due 
to star formation (Ehlerova and Palous, 1996). In a new search by Ehlerova 
and Palous (2004), more than 600 shells have been identified in the Leiden- 
Dwingeloo HI survey of the Milky Way. In Fig. E|we show the re-identification 
of a shell previously discovered by Heiles (1979) and a newly discovered shell. 
The distribution of them in the radial galactocentric direction and in the 
direction perpendicular to the galactic disk is similar to stellar distribution 
supporting the idea of a connection between massive stars and shells. 

11.2 The collapse of shells 

We discuss the supersonic expansion of shells and sheets from regions of local- 
ized deposition of energy and address the question if and when they fragment 
and collapse due to gravitational instability. The energy input from an OB 
association creates a blast-wave which propagates into the ambient medium 
(Ostriker & McKee, 1988; Bisnovatyi-Kogan & Silich, 1995). The schematic 
representation of the situation is shown in Fig. After the initial fast ex- 
pansion the mass accumulated in the shell cools and collapses to a thin struc- 
ture, which is approximated as infinitesimally thin layer surrounding the hot 
medium inside. Neglecting the external pressure and assuming the constant 
energy input L, the self-similar solution for radius R, expansion velocity V 
and column density S s /, is (Castor et al. 1975; Ehlerova & Palous, 2002): 

=53-1 x ( 1Q51 crg L Myr -^ x x (^)*pc (4) 

V (t) =31.2x ( 1Q51 crg L Myr -^ x x (^j^kms- 1 (5) 
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Figure 3: Re-identification of the shell GS 128+01-105 discovered by Heiles 
(1979) - upper frame, and a newly discovered shell - lower frame. 
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Figure 4: Schematic representation of a supershell expanding around an OB 
association. 



1 3 

£(t)* = 0.564 x( 1Q51 cr L eMyi .^ y x (^- 3 ^y x (j^Y M pc- 2 , (6) 

where n is the density, /j, is the mean atomic weight of the ambient medium 
and t is the time since the beginning of an expansion. 

The linear analysis of hydrodynamical equations including perturbations 
on the surface of the shells has been performed by Elmegreen (1994) and 
Wiinsch & Palous (2001). The fastest growing mode is: 

where c a h is the sound speed inside of the expanding shell. 

In Fig. we give the time evolution of the fastest mode. At early times, 
for t < U, the shell is stable. t& is the time, when the fastest mode starts to 
be unstable: 

t b = 28.8 x ( -2*^) 1 x ( —± -\ " x ( JL " ) "* My, (8) 

Vkms" 1 / VlO 51 erg Myr" 1 / V1.3cm- 3 / J w 

Later, for t > tb, when the expansion slows down and reduces the stretching, 
which acts against gravity, and when the shell column density increases, the 
shell starts to be gravitationally unstable. For ambient densities similar to 
values in the solar vicinity, n ~ 10 _1 — 10 2 , tb is a few 10 7 yr, which means 
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Figure 5: The fastest growing mode. 



that the gravitational instability is rather slow compared to turbulent collision 
times and galactic differential rotation, tb is much smaller in high density 
medium of GMC and dense cores, where it is ~ 10 4 yr only. Thus the shell 
gravitational instability is particularly important inside the GMCs. 
The dispersion relation of the shell gravitational instability is: 



where 77 is the dimensionless wavenumber and A is the wavelength of the 
perturbation: 77 = 2nR/X. It is shown in Fig. HO it gives the wavelength 
interval of unstable perturbations. 

The resulting number of fragments is inversely proportional to the frag- 
ment growth time t growth — Jf^j) ■ Rapidly growing fragments are more 
frequent in the final mass spectrum than the slowly growing fragments. 



N = lo /x //n3 . (10) 



Thus the number of fragments in a given volume of radius R is 

_R 3 

(A/4) 

A fragment with the wavelength A has the mass 

m = ^tt(A/4)V (11) 

We derive the mass spectrum £(m) = j^: 

4 

£,(m) = n R 3 p uj m- 2 . (12) 
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Figure 6: The dispersion relation 



If the dispersion relation ui(rj) were a constant, the slope of the mass 
spectrum £(m) oc mT a would be exactly approximated with a power law 
slope a= 2. But in our case u> is not only the function of rj but also of the 
time t. 

We assume that t growt h for a given r\ is inversely proportional to the aver- 
age value of oj for this r\ since the time %{rj) when a given mode starts to be 
unstable. The time average Q for a given rj is calculated using the equation: 



The resulting mass spectra for different values of n, c s h and L, as they have 
been derived using the thin-shell approximation (4) - (6), are shown in Fig. 
The high mass parts are well approximated by the power law with a slope 
a = 2.2 -2.4. 

In low density medium, the collapse of expanding shells forms fragments 
with masses comparable to GMC, however, the collapse time is rather long, a 
few 10 7 yr. In high density medium of GMC cores, the collapse time is rather 
short, ~ 10 4 yr, the masses of individual fragments are close to stellar and the 
initial mass function of fragments has a power law slope close to the Salpeter 
(1955) value -2.35. We conclude that the fragmentation of expanding shells 
qualifies as a possible process triggering the star formation in environments, 
where the density is high enough, or where it has been increased due to mass 
accumulation. 
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Figure 7: The mass spectrum of fragments of an expanding shell at time 5ib- 

12 Formation of super-star clusters 

Young and massive (10 5 — 10 7 M©) super-star clusters are observed along 
tidal arms and bridges of the colliding galaxies. A model of a star form- 
ing factory is proposed by Tenorio-Tagle et al. (2003). This model invokes 
pressure-bounded, self-gravitating, isothermal cloud (Ebert, 1955; Bonner, 
1956), which becomes gravitationally unstable when sufficiently compressed. 
The gravitational instability allows the cloud to enter isothermal collapse. 
As the collapse proceeds a first generation of stars is formed in the center 
of the cloud. The mass and energy feedback of the first generation of stars 
has an important impact on the collapsing cloud. Stellar winds and super- 
novae compress the infalling material forming a dense shell (see Fig. |HJ. The 
shell is able to trap the ionizing radiation and winds of the first generation 
of young stars. At the same time the shell fragments forming stars with a 
high efficiency. Thus the cloud, which is compressed in the case of the galaxy 
encounter, forms a new super-star cluster. 

13 Evolution of the star formation rate 

The investigation of colors in Hubble deep fields (Madau et al., 1996; 1998; 
Rowan-Robinson, 2003) provides extinction as a function of z: it has been 
higher at z = 0.5 — 1.5 than locally, and lower at z > 2. Related models of the 
star formation history show the steep decline of the star formation rate since 
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Figure 8: The shell inside of a collapsing cloud 
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z = 1 , the present star formation has in average about an order of magnitude 
smaller rate compared to the level at z = 1 — 1.5. Even deeper in the past the 
correction for the dust extinction remains highly uncertain and a conclusion 
on the evolution of the star formation rate for z = 2 — 6 is difficult. 

This opens the question what drives the star formation. Bars in the central 
parts of galaxies does not seem to change substantially between z = 1 and 
present. They do trigger the star formation, however, there has to be more 
partners in the game. The fuel - the gas - is depleted in galaxies not only due 
to star formation, but also since the ISM is removed from the galaxy disks 
with tides, with the ICM ram pressure stripping and with the star formation 
feedback. Number of galaxy interactions also decreases. We conclude that 
the decline of the star formation is due to lower rate of triggering from galaxy 
interactions in a combination with the starvation, since the amount of gas 
in the star forming disks is reduced by environmental effects and by the star 
formation itself. 
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